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Abstract. The numerical analysis of strongly interacting nanostructures requires powerful techniques. Re- 
' cently developed methods, such as the time-dependent density matrix renormalization group (tDMRG) 

approach or the embedded-cluster approximation (ECA), rely on the numerical solution of clusters of fi- 

■ nite size. For the interpretation of numerical results, it is therefore crucial to understand finite-size effects 
in detail. In this work, we present a careful finite-size analysis for the examples of one quantum dot, as 

, well as three serially connected quantum dots. Depending on "odd-even" effects, physically quite different 

results may emerge from clusters that do not differ much in their size. We provide a solution to a recent 
controversy over results obtained with ECA for three quantum dots. In particular, using the optimum 
clusters discussed in this paper, the parameter range in which ECA can reliably be applied is increased, 
as we show for the case of three quantum dots. As a practical procedure, we propose that a comparison of 

' results for static quantities against those of quasi-exact methods, such as the ground-state density matrix 

renormalization group (DMRG) method or exact diagonalization, serves to identify the optimum cluster 

■ type. In the examples studied here, we find that to observe signatures of the Kondo effect in finite systems, 
the best clusters involving dots and leads must have a total ^-component of the spin equal to zero. 

PACS. 73.63.-b Electronic transport in nanoscale materials and structures - 73.63Kv Quantum dots - 
71.27.-|-a Strongly correlated electron systems 



1 Introduction 

The prediction and experimental observation of the Kondo 
effect in quantum dots [1-4] and single-molecule conduc- 
tors has stimulated considerable interest in strongly cor- 
related nano-scale systems, as discussed in recent reviews 
[5-7] . The related experimental and theoretical efforts are 
not only motivated by novel emerging physical phenom- 
ena, but also by the possible technological applications of 
nanodevices, and in particular by their transport prop- 
erties. Experimental results for the conductance through 
nanostructures have been reported for single dots, side 
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and linearly coupled two dots, as well as small molecules 
[1,8,9]. 

Stimulated by the rich physics harbored by interact- 
ing nanostructures, the field has seen a rapid develop- 
ment of powerful numerical techniques to obtain the con- 
ductance. These include the well-established numerical- 
renormalization group approach (NRG) [10-13], various 
density-matrix renormalization group (DAIRG) based meth- 
ods [14-16], including in particular, the recently devel- 
oped time-dependent DMRG (tDMRG) [17-22], Quan- 
tum Monte Carlo simulations (see, e.g., Ref. [23]), flow- 
equation approaches [24] , as well as exact diagonalization 
combined with an embedding procedure, the embedded- 
cluster approximation (ECA) [25-28]. We wish to distin- 
guish the latter method from other cluster-embedding ap- 
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preaches used in this field, sometimes dubbed the embed- 
ding method [29-31]. 

Beyond numerical methods, there are many analytical 
techniques to the problem of transport in quantum dots, 
ranging from perturbative schemes (see, e.g., Refs. [32- 
34]), real-time renormalization methods [35], and func- 
tional renormalization group approaches [36,37], to slave- 
boson calculations [38], and variational methods [39,40], 
just to mention a few. Furthermore, density functional 
theory (DFT) is widely used to model molecular conduc- 
tors (see, e.g., Ref. [41] for a review). Such DFT-based 
approaches successfully account for the coupling of the 
leads to extended molecules by tackling the molecule plus 
a substantial part of the leads (see, e.g., Refs. [42,43]). In 
that sense, there is a similarity with the embedded-cluster 
approximation, as will become more obvious in the dis- 
cussion of EGA in Sec 2.2. While DFT based methods 
have been shown to correctly describe Coulomb-blockade 
physics [42], interacting resonant level models [44], and 
finite-bias phenomena [43,45], the effects of strong corre- 
lations as manifest in the Kondo effect still seem to be 
out of reach for current DFT schemes. In that sense, DFT 
based methods are complimentary to other approaches 
mentioned here that emphasize the effect of strong cor- 
relations. 

Time-dependent DMRG and EGA are particularly pro- 
mising methods to address complex nanoscale systems 
such as molecular conductors. Both methods can in prin- 
ciple include the modeling of phonons, and we mention, as 
an example, the EGA study Ref. [46] on a molecular con- 
ductor with a center-of-mass motion. Other applications of 
EGA on timely problems of experimental significance in- 
clude studies of nonlocal spin control in quantum dots [47] , 
the SU(4) Kondo effect [48], and multilevel quantum dots 
[49]. Furthermore, it is one of the advantages of DMRG 
[50-53] and EGA that spatial correlations can be analyzed 
not only within the interacting region, but also throughout 
the entire system. Therefore, it is possible to study how 
the screening of the magnetic moment of quantum dots 
occurs, which relates to the notion of the so-called Kondo 
cloud [54-57]. The computation of extended spatial corre- 
lations within NRG has become possible only recently [58] . 
Moreover, existing codes for both EGA and tDMRG can 
easily be adapted to different problems, and EGA calcu- 
lations can be performed at low computational costs, i.e., 
the requirements in GPU time needed are rather mod- 
est compared to DMRG and NRG. This is, compared to 
NRG, of particular advantage in the case of multi-channel 
problems, such as realized in the case of three quantum 
dots. 

As both DMRG and EGA are based on real-space 
schemes, for the interpretation of numerical results, it is 
crucial to understand how subtle many-body effects, such 
as the Kondo effect, manifest themselves on finite sys- 
tems. It is the purpose of this work to present a detailed 
analysis of the properties of strongly interacting nanos- 
tructures and their finite-size scaling, in both DMRG and 
EGA calculations. Our work first provides a discussion of 



(a) One quantum dot 




(b) Three Q-dots 
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Fig. 1. Sketch of the models studied in this work: (a) a single 
quantum dot embedded into noninteracting leads; (b) three 
serially coupled dots embedded into leads. Full symbols denote 
dots (namely, sites with a finite onsite Coulomb repulsion U > 
0), while open symbols represent tight-binding sites. 



the one-dot case and then covers three serially coupled 
quantum dots as well. The models are depicted in Fig. 1. 

We discuss how the emergence of Kondo physics is 
reflected in static properties, such as spin and charge fluc- 
tuations and spin-spin correlations. In addition, we dis- 
cuss the coirductance of these structures, obtained with 
the EGA method, as a function of model parameters and 
gate voltage, which controls the dots' filling. As our main 
result, we argue that it is important to consider the ge- 
ometric properties and global quantum numbers of finite 
systems in order to correctly interpret the results from ex- 
act diagonalization (ED)-based approaches such as EGA. 
The main conclusions are relevant as well for real-time 
simulations within tDMRG performed with tight-binding 
leads. We mainly gauge our results against the most es- 
tablished methods, the Bethe ansatz solution of the single 
impurity problem [59,60], sum rules valid for Fermi liquids 
[61], and, for multi-dot structures, NRG [13]. 

The existence of odd-even effects in NRG calculations, 
which use logarithmically discretized leads (sometimes 
caUed Wilson chains, see, e.g., Refs. [16,22]), has been 
emphasized by Wilson [10], with the important observa- 
tion that the renormalization procedure of NRG on Wil- 
son chains results in two different fixed points, depending 
on whether the number of NRG iterations is odd or even. 
In this work, we discuss odd-even effects that emerge in 
the finite-size scaling of, e.g., the conductance, when the 
dot is coupled to leads in a real-space representation [62] . 
Such even-odd effects and, more generally, the relevance of 
boundary conditions for transport, have been studied be- 
fore, e.g., for small chains of dots coupled to wires [63], or 
quantum box problems [64] . The purpose of our work is to 
discuss these effects with a perspective on the two methods 
of interest here, i.e., tDMRG and EGA, and most impor- 
tantly, on the calculation of conductances for the single- 
impurity Anderson impurity model and a three-quantum 
dot model. We mention that even-odd effects and their 
influence on the conductance have previously been stud- 
ied by static DMRG in Ref. [29] . The model studied there 
was spinless fermions interacting in a constrained region 
embedded in noninteracting leads in a ring geometry. 

The case of three dots [27,28,36,40,65-70] has attracted 
considerable attention from the theoretical side. Most stud- 
ies find perfect conductance through the dots at half fill- 
ing, independently of model parameters such as hopping 
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matrix elements between the dots and onsite Coulomb re- 
pulsion [36,40,65-69]. This result is at odds with earlier 
results obtained from EGA since a conductance dip has 
been reported in Ref. [28]. Here, we resolve this issue in 
favor of the picture promoted in Refs. [36,40,65-69] and 
show that the controversial dip is not due to a general 
incapability of EGA to capture the essential physics, but 
due to quite subtle finite-size effects which have not been 
appreciated in some earlier EGA calculations [28]. 

The EGA method incorporates the numerical solution 
of a small cluster that contains the interacting region - 
e.g., a quantum dot - by means of exact diagonalization. 
One example of a relevant global quantum number of such 
small clusters is S^^^^^, the z-component of the total spin. 
We illustrate that in certain cases, clusters with an overall 
•^totai 7^ hsdi filling, and with open boundary condi- 
tions may exhibit qualitatively different properties from 
those with S^^^^^ = 0. Furthermore, as a practical proce- 
dure, we propose that a quantitative comparison of static 
properties, such as charge variations with gate potential, 
between EGA and ground-state DMRG or the Lanczos al- 
gorithm, allows us to determine the optimum cluster-type 
for a given model. A Lanczos solver is part of standard 
EGA codes as we shall explain later. 

The Hamiltonian common to the systems analyzed here 
consists of three parts: the noninteracting leads -fficads, 
the coupling between the interacting region and the leads, 
iJhy, and the interacting region described by -ffint: 

H = i?lcads + ^^hy + -f^int , (1) 

where the leads are described by: 

Hlcads = X! \-4,a'^W.cr + h.c] 

1=1 

't E [cLc/+i.a + h.c] . (2) 
;=Wi,+JVi„t+i 

As usual, c[^l denotes a fermion annihilation (creation) 
operator acting on site I, with a spin index a = |,J,. 
Summation over a repeated index a is implied throughout 
the paper, ni^^ = c\ ^ci^„ is the local particle operator 
with spin a and ni = n/j -|- ni^i is the particle number 
(or charge) operator on site I. The total system size is 
TV = Nl + A^int + A'fi, where Ni^i^j^-^ is the number of sites 
in the left(right) lead and A'int is the number of interacting 
sites in the center of the system. The hybridization term 
can be written as: 

^hy = -t'^ic'^^^^CNr^ + ha + cl ,^C^+l,„ +h.C.) , (3) 
a 

with X = Nl + Nint- Unless otherwise stated, i = 1 is 
the unit of energy. The quantum of conductance in this 
notation is thus Gq = 2, where the factor of two is due to 
the two spin channels. 

The rest of the paper is organized as follows. In Sec. 2, 
we introduce the two methods used here, DMRG and 



EGA. In Sec. 3, we revisit the case of a single quantum dot. 
We point out that tDMRG calculations result in quite dif- 
ferent behaviors, depending on the type of open clusters. 
We qualitatively discuss spin-spin correlations and calcu- 
late the conductance as a function of gate voltage with 
EGA. In Sec. 4, we present our numerical results for static 
properties of three quantum dots obtained from DMRG 
and EGA calculations, with a special focus on their finite- 
size dependence. This allows us to identify the optimum 
clusters for the calculation of dynamic properties and the 
conductance. Section 5 provides a summary and conclu- 
sions. 

2 Methods 

2.1 Time-dependent DMRG 

For ground-state DMRG calculations, we use the finite- 
system size algorithm [52] and we give the number of 
states M used in either the figures' captions or in the 
text. The set-up for our adaptive time-dependent DMRG 
[71,72] calculations of transport properties has been de- 
tailed in Ref. [17]. Here we just repeat that the current 
J(t) [where r denotes time in units of 1/t] is measured 
as the current across the link between the leads and the 
interacting region [see Eq. (3)]: 

jx,x+l = -it'J2^'^x+l,aCx,a " h.C.) 

J = (0'jVt+JVi„t,JVt+Ari„t+i)T + {jNL,NL + l)T)/'2 (4) 

with X = Nl + Nint or X ~ Nl- ( • )t denotes the expecta- 
tion value taken in the time-dependent state. The external 
bias - an onsite potential zLAVui applied to the leads - 
is AV ~ 10~^t. We use a Trotter-Suzuki break-up of the 
time evolution operator in our adaptive tDMRG scheme 
with a time step of 6t ~ 0.05 [71,72]. The discarded 
weight during the time evolution is kept below 10"^*^. For 
the tDMRG data to be discussed in Sec. 3.1, M < 800 
states needed to be kept to meet this requirement. With 
this choice of tDMRG parameters, the real-time data ob- 
tained with a discarded weight of 10^^° and 5 • 10~^° turn 
out to be indistinguishable from each other on the scale 
of the plots. Open boundary conditions are imposed in all 
DMRG calculations. As for the notation adopted through- 
out the manuscript, for our DMRG results, ( • ) denotes 
the expectation value taken in the ground-state calculated 
with the method in a subspace given by a fixed total num- 
ber of electrons and a fixed value of S^^^.^^. Further note 
that on chains with an odd number of sites and at half fill- 
ing, the results on static properties and real-time currents 
to be presented in the following sections do not depend on 
whether we work at S'^jai — or Sfotai = ^1/2, with 
the trivial exception of (Sf), which changes its sign. For 
N odd, we hence work in the S^^^^^ = 1/2 subspace only. 

2.2 EGA 

The EGA method relies on the numerical determination 
of the ground-state of a cluster with open boundary con- 
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ditions. In the following, we briefly sketch details of the 
method. 

The EGA method tackles the impurity problem in three 
steps. First, the infinite system is artificially cut into two 
parts: one part C that contains the interacting region 
plus as many noninteracting sites of the leads as possi- 
ble (this part will, from now on, be referred to as 'the 
cluster'), and a second part R (the 'rest'), consisting of 
semi-infinite chains positioned at the left and the right in 
relation to the cluster C. The number of sites in C is de- 
noted by -/Ved- Second, Green's functions (GF) for both 
parts are computed independently: current implementa- 
tions of EGA utilize the Lanczos method to calculate the 
interacting region's GF, while those of the part R, being 
noninteracting, can be computed exactly as well. In a final 
step, the artificially disconnected parts are reconnected by 
means of a Dyson equation, which dresses the interacting 
region's GF. This step, the actual embedding, is crucial 
for capturing the many-body physics associated with the 
Kondo effect. Moreover, although the clusters that can be 
solved exactly by means of a Lanczos routine are rather 
small, being of the order of A^ed ~ 12 sites only, the em- 
bedding step largely compensates for that by dressing the 
cluster GF and effectively extending the many-body cor- 
relations, induced by the presence of the impurity, into the 
semi-infinite chains R. 

Wc now further detail these steps. The Hamiltonians of 
the left and right semi-infinite, tight-binding chains, i.e., 
the noninteracting R part, are described by 



with 



oc 

Hsc-R = -t (cj'^q+iCT +h-c.), (5) 

1=Neu + 1-X' 

where in this notation, the sites labeled by i = 1, . . . , Aed 
are inside the cluster C. The semi-infinite chains are con- 
nected to the cluster by the following term: 



(6) 



where V = t is the hopping in the broken link, connect- 
ing parts R and C, used for the embedding procedure. 
The GF for the cluster C and for the semi-infinite chains 
are calculated at zero temperature. Fixing the number of 
particles m and the z-axis projection of the total spin, 
"^totai' ^^"^ ground state and the one-body propagators be- 
tween all the clusters' sites are calculated. For example. 



9'. 



('".•S'total) 



the undressed Green function for the cluster. 



propagates a particle between sites i and j inside the clus- 
ter. For the noninteracting, semi-infinite chains, the GFs 

and gN^n+i ^^^'■^ sites and A^ed + 1, located at the 
extreme ends of the semi-infinite chains, at left and right 
to the cluster, can be easily calculated as well. 

The Dyson equation to calculate the dressed GF ma- 
trix elements G,-™''^'"*"'^ can therefore be written as 



G 



(■'"'■Stotal) 



9i,j 



„(™.'S'to, 

9i.i 



„(™''St"otal) T/ /^('"'■St'otal) 

~yi,NEn ^ '^Neu + IJ 



(7) 



^'('"■■^total) 
r'(™.'St'otal) 



V Gi' 



9^ 



9Ned + 1 



('".■Stotal) 



T/ /^(™.'S't~otal) 

^ ^NEn,3 



(8a) 
(8b) 



where V, as mentioned above, is defined according to 
iJicst, Eq. (6). Equations (7) and (8) correspond to a chain 
approximation, where a locator-propagator diagrammatic 
expansion is used [73,74]. Note that EGA is exact in the 
case of [/ = 0. 

As mentioned before, the calculation of the propaga- 
requires that fixed quantum numbers m and 



tor gl ^ 



'5'totai tie used. However, after the embedding procedure, 
these quantum numbers are not good quantum numbers 
for the cluster anymore. Therefore, we have to incorpo- 
rate processes into the EGA method that allow for charge 
fluctuations in the cluster C. To accommodate this re- 
quirement, different implementations of EGA have been 
devised, either by including different spin-mixing strate- 
gies [25,47,75,76] or by moving the Fermi energy of the 
leads [77,78]. 

For this work, we will adopt a spin-mixing strategy dif- 
ferent from the one used in a previous work analyzing the 
three-dots geometry [28] . We will show that this different 
strategy is instrumental in decreasing finite-size effects ob- 
served in Ref. [28]. The spin mixing proceeds as follows. 
First, a cluster GF with mixed charge is defined through 



9. 



(™+P.P'S'total) _ 



■1,3 



= (1-P) 9l 



(?T1,0) 

3 



P5ir'''^'"'' (9) 



where p takes values between and 1, and we assume that 
m is even, in which case, the corresponding S^^^^i = 0. 
In addition, note that for the cluster with charge m + 1, 
'5'totai takes values ±1/2. The matrix element g^^J^P'^^tot^i) 
corresponds to a situation where the charge in the cluster 
is between m and m + 1. The total charge in the cluster, 
before embedding, can be easily calculated as 



(p) = (1 ^ P)™ + P(™ + 1) ~ m + p. (10) 



Using Eqs. (7) and (8), the dressed GF G("'+P^pSt'otai) is 
obtained, and from this result, the charge in the cluster 
can be calculated: 



QpSL...(^p) = Jl r"lm{^G 



{uj)}duj, 



(11) 

where Ep is the Fermi level. The value of p is calculated 
self-consistently, satisfying 



(12) 



If there is spin reversal symmetry, e.g., no magnetic field 
is applied, one can calculate the total GF as 



G 



(P)-1 E G,5+^'^^— \ (13) 

Sfotal = ±l/2 



where p satisfies Eq. (12). It is important to emphasize 
that the charge fluctuations taken into account by Eq. (9) 
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arc the ones between the cluster and the rest of the system, 
and not just the ones at the interacting region described by 
iJint • The latter ones involve a very localized neighborhood 
of the dot and as a consequence, are typically already well 
described on isolated clusters only. 

As will be shown in Sec. 4.2, this alternative way of 
obtaining the dressed GF leads to improved results in the 
case of three dots. Note that in the case of one dot, no 
substantial difference in the final results obtained from ei- 
ther the old spin-mixing strategy [25,47,75,76] or the one 
described here is found. Yet, for the sake of consistency, 
wc apply the same way of spin mixing as described here 
throughout the paper. 

It is noteworthy to point out that the self-consistent 
solution for the charge-mixing parameter p is either or 
1 in the Kondo regime. Therefore, in the Kondo regime, 
where charge fluctuations are suppressed, no charge mix- 
ing is done at all. The parameter p takes a finite value 
< p < 1 as the gate potential drives the system into the 
mixed-valence regime, and the purpose of the charge mix- 
ing is thus mainly to smoothen out the transition from an 
m electron to an m ± 1 electron ground state, that for the 
bare cluster would be a discontinuous one. From a con- 
ceptual point of view, the introduction of the parameter 
p establishes a necessary consistency in the theory guar- 
anteeing that the number of electrons in the cluster is the 
same before and after the embedding process. 

For completeness, we give the actual expression that 
we compute the conductance G from [25,26,79]: 

G = Go [t^Plc.ds{EF)?\G,,r{EF)\''. (14) 

Here, picadsiEp) is the density of states of the leads at the 
Fermi-level Ep and Gi^r{'-^) is the dressed Green function 
that propagates an electron from the first site / left of the 
cluster to the first site right r of the cluster. 



3 Transport through a single quantum dot 
revisited 

In this section, we focus on the application of tDMRG and 
EGA to the study of the conductance of a single quantum 
dot. The interacting portion of the Hamiltonian (1) is: 

i^int = t^"dot,T'T'dot,i + ^9"'dot , (15) 

where the Goulomb repulsion U represents the charging 
energy, Vg is the gate potential, and ndot = n^ot,-] +"-dot,i 
is the dot's charge. The model is particle-hole symmetric 
at Vg = -C//2. 

As has been shown in Ref. [17], tDMRG is capable 
of producing zero-temperature transport properties of the 
single-impurity Anderson model for U/{2t'^) < 4, such as 
the conductance as a function of gate potential or its de- 
pendence with a magnetic field. In that previous work, 
it was noticed that the behavior of real-time currents on 
finite systems with open boundary conditions is quite dif- 
ferent depending on odd-even effects. When the number 
of sites in each lead is even or odd, we have a S^^^^y = 1/2 




5 10 15 20 



time X 

Fig. 2. One dot, t' /t = 0.4, Vg = -17/2. Main panel; current 
J{t)/AV vs. time for = 31 (oio), A = 32 {ole), and N = 
33 {ele) at U /t = 2 (tDMRG). Inset: current J{t)/AV vs. 
time at [/ = for N = 31,32,33 (exact diagonalization) . We 
divide J(t) by the bias {AV = lQ~^t). Perfect conductance 
corresponds to J{t)/AV = 2 (indicated by horizontal, dotted 
lines). 

ground-state, while with an odd number of tight-binding 
sites in one lead and an even number in the other, the 
system has an overall S^^^^^y = ground state. Here we 
present a detailed scaling analysis of static properties that 
explains why the latter type of clusters is preferable, as 
suggested in Ref. [17]. 

Finally, we discuss finite-size effects in the context of 
EGA, which are of similar origin, and we provide a quan- 
titative comparison with exact results obtained through 
the Friedel sum rule [61] and ground-state DMRG. 



3.1 Motivation: tDMRG results 

Let us motivate our study by discussing results from tDMRG 
for the current through a quantum dot at small external 
biases. On an open chain with one embedded dot, we an- 
alyze three configurations: (i) an odd number of sites in 
both leads (oJo); (ii) an even number of sites in both leads 
{ele), (iii) or an overall even number of sites {ole). We re- 
peat some of the main results obtained in Ref. [17] for the 
ole cluster-type, which has mostly been used in that work. 
First, both at zero and finite U , the current is almost con- 
stant in time over an extended period in time, allowing 
us to assign a conductance value G{N) = J{t)/AV to 
each chain of finite length. Note that the length in time of 
such constant current is limited by twice the time that it 
takes the current to be reflected at the open boundaries. 
Second, performing a finite-size scaling of G{N), one finds 
that it extrapolates, following a simple scaling law of the 
form G{N) ~ 1/N, to the known exact results from Fermi 
liquid relations [61] for U/2t''^ < 4. Third, ele clusters 
have been found to exhibit a slower convergence and gen- 
erally lower current values. Our discussion in this section 
will mostly focus on clusters with S^^^^^ = 1/2, i.e., ele 
and olo ones. 
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The behavior of a current driven by a small bias is 
illustrated in Fig. 2 for clusters of an intermediate length 
and at half filling. It is important to notice that we choose 
values for U and t' as well as cluster sizes such that the 
finite-size effect we wish to address is clearly visible. As a 
drawback, obviously, a steady state with G ~ Go is not 
reached for the three clusters and system sizes considered 
in Fig. 2. Longer chains are needed for this set of U and 
t' since (i) finite-level spacings possibly cut off the Kondo 
resonance and (ii) the time-scale for reaching the steady 
state [80,81] may exceed the time at which the current 
reverses its sign. 

We now turn to this section's main point, the discus- 
sion of even-odd effects in the current, ele clusters result 
in the smallest current, ole clusters yield a larger current 
and exhibit the aforementioned behavior of being approx- 
imately constant in time ( J(t) «const) before the reversal 
of the current's sign occurs, olo clusters allow for a sub- 
stantially larger current that even exceeds - if considering 
the conductance G = J/AV - perfect conductance Go = 2, 
but that is never constant in time. These differences are 
quite striking, considering that the system sizes used in 
Fig. 2 differ by only one site in the leads. It is worth point- 
ing out that these strong finite-size effects arc in principle 
not caused by the presence of the Coulomb repulsion on 
the dot, but can even be seen in the noninteracting case 
U = 0, as displayed in the inset of Fig. 2. Note that when 
the system size is increased, J{t)/AV measured on olo 
clusters approaches perfect conductance from above. We 
stress that the strong differences in J(t) depending on the 
cluster type are not due to current-ringing [80,81], which 
is most pronounced away from half filling {Vg ^ —U/2) or 
at large biases [18,71,82]. 

The crucial difference between ole on the one hand 
and the other two clusters on the other hand is the finite 
magnetic moment S^^^^^^ = 1/2 present in the system in the 
case of an overall odd number of sites. As we shall see, this 
finite moment is not homogeneously distributed over the 
system. Let us first explain the suppressed conductance 
seen in the case of the ele clusters. Figure 3(a) shows the 
local spin density (S^) for N = 17 sites and several values 
oiU/t ^ 0, 1, 2, 4. Clearly, for any U, most of the S^^^^i is 
carried by the central site, i.e., the quantum dot [62]. The 
reason for this inhomogeneous distribution is easy to un- 
derstand, considering the critical tendency of ID systems 
towards antiferromagnetic (AFM) correlations at half fill- 
ing. For that reason, the spin cannot be homogeneously 
distributed over all sites. We further observe that on open 
systems, the spins on the first and last tight-binding site 
point up, which then, with an even number of sites to the 
left and right of the dot, favors the spin on the dot to 
point up as well. In other words, both leads form an ap- 
proximate singlet, leaving the remaining moment S^^^_^l on 
the dot. This is strictly true in the limit oi t' = 0, where 
(■^dot) = 1/2, while each lead has a vanishing (S'^tai)- 
creasing t' then reduces (S^^^). Upon increasing [/ at a 
fixed t' > 0, a larger fraction of S^^^^^ is pinned at the dot 
site, as charge fluctuations that tend to reduce the local 
moment are suppressed by increasing U. 




5 10 

site i 



Fig. 3. One dot, t' /t = 0.4, Vg = -U/2. Spin density (Sf), 
for (a) A'^ = 17 sites (ele) and (b) TV = 15 sites (olo) with 
U/t = 0,1,2,4. (a): Even at [/ = 0, the z-coniponent of the 
total spin is essentially carried by the central site, i.e., the 
dot. (b): In the olo case, overall, every second site carries a 
positive (Si). The dotted vertical lines denote the position of 
the dot. DMRG data with M = 300. The results presented 
here for olo and ele agree with Ref. [62]. DMRG calculations 
are all performed in subspaces with a fixed number of fermions 
and Sfotai, and the symbol {■), here and throughout, refers 
to taking the expectation value in the ground state of such 
subspaces. 

Qualitatively, the situation of an ele cluster is equiva- 
lent to a quantum dot in the presence of a magnetic field: 
locally, the Z2 symmetry S^^^ <-^- — S'^^t is broken on ele 
clusters. A magnetic field is well known to split the Kondo 
resonance thus leading to a reduced conductance [61]. In 
order to support this picture, we use ele clusters and ap- 
ply a real magnetic field onto the dot: 

-fffiold = /jdotS'dot; 'S'dot = 2('^dot,T - '^dota) • (16) 

As expected, when {S^^^) decreases upon applying the 
field h^ot ~ as now the field favors a negative spin in the 
dot - the current increases, as we show in Fig. 4. For the 
parameters of Fig. 4 {N = 33, ele) we infer from the fig- 
ure's inset that the spin projection on the dot vanishes at 
about /idot ~ O.SAt. Consistently, the current measured in 
tDMRG simulations at h^^t ~ 0.35t reaches values similar 
to those observed using olo clusters of comparable system 
size. This establishes a connection between the dot's po- 
larization and the achievable currents. Note further that 
ele clusters behave exactly like olo ones if the filling is 
kept at (iV± l)/iV (not shown). 

Finally, note that the finite spin projection in the clus- 
ter's ground-state is an issue of parameters. Increasing the 
coupling to the leads, i.e, t' , at fixed U and system size 
(say, N = 33, ele), delocalizes the spin and a larger cur- 
rent is then measured in tDMRG calculations. 

For the second type of clusters with a nonzero S^^^^^^, 
olo, we find that the spin is equally distributed over every 
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time T 



Fig. 4. One dot, t' jt = 0.4, U jt = 2, Vg = -C//2. Main paneh 
tDMRG data for the current J(t)IAV vs. time for iV = 33 
and several values for a magnetic field /idot applied to the dot: 
/idot/t = 0, 0.1, 0.2, 0.35 [see Eq.(16)]. The inset shows (Sdot) as 
a function of /idot for A'' = 33 and 65, calculated with ground- 
state DMRG (M = 300). 

second site [62], displayed in Fig. 3(b). Since the out-most 
sites carry a positive {Sf) > 0, the dot now has a nearly 
vanishing {Sl^^). Consistent with the notion of {S^^^) w 0, 
the spin projections on the dot and all other sites are 
practically independent of the Coulomb repulsion U. 

It is important to stress that at fixed odd N, both 
types of behavior can be produced: when shifting the dot 
in a = 33 cluster by one site, a 17-dot-15 configuration 
is obtained. The current now behaves almost exactly like 
the one measured on an TV = 31 {olo, 15-dot-15) clus- 
ter and so does the z-component of the spin. The same 
holds for an asymmetric = 31 cluster with a 16- dot- 14 
configuration: the dot now carries most of the S^^^^^ and 
the current is suppressed, similar to the case of A'^ = 33 
{16-dot-16, ele). The corresponding DMRG results are 
not shown in the figures, but are quantitatively similar to 
Figs. 2 and 3. These considerations show that indeed the 
distribution of S'tot^i i which depends on the geometry and 
model parameters, affects the conductance. 

We finally emphasize that the issue of a non-uniformly 
distributed {S^) on half- filled chains with an odd num- 
ber of sites and constant hopping matrix elements in the 
leads also emerges in other interacting systems that do not 
feature Kondo physics, such as transport through a Mott 
insulator or an even number of quantum dots [83] . Its in- 
fluence on the finite-size scaling of transport properties is 
most significant in Kondo problems, though. 

3.2 Static properties: Gate potential dependence 

To render the line of reasoning outlined in Sec. 3.1 more 
quantitative we proceed with a systematic discussion of 
static properties that are crucial to characterize a dot in 
the Kondo regime: spin and charge fluctuations on the 
dot, as well as spin-spin correlations between the dot and 



_ (a) One dot, t7t=0.4, U=4t 




Fig. 5. One dot, static properties computed with DMRG as a 
function of gate potential Vg/U at t' /t = 0.4, U/t = 4, compar- 
ing different cluster types: olo {N = 15, circles), ole {N — 16, 
squares), and ele {N =17, stars). Triangles in (a) and (b) are 
for a large system of TV = 200 sites (M = 300 DMRG states), 
(a) Charge fluctuations (5ndot)^; (b) spin fluctuations (dS^^^)^; 
(c) local spin density {S^ot) on the dot; and (d) spin-spin cor- 
relations S^dot between the dot and the first neighboring site 
in the leads. For ole clusters such as A'^ = 16, there is an asym- 
metry between the left and the right lead: the dot is mostly 
screened by the left lead, i.e., the one that has an odd number 
of sites. 



the leads. For this purpose, we focus on the parameter set 
t'/t = 0.4 and U/t = A and the dependence on the gate 
potential, keeping the full system at half filling. 

First of all, let us recall some of the hallmark features 
of the single-impurity Anderson model at T = [13,61]: 
(i) a slowly varying charge (ndot) ~ 1 on the dot in a 
broad gate-potential window — 1 < Vg/U < for T <g; [/, 
where F denotes the hybridization parameter. This prop- 
erty, i.e., the almost constant charge, becomes more pro- 
nounced as U I r increases, (ii) Enhanced spin fluctuations 
on the dot. (iii) Large charge fluctuations at the charge 
degeneracy points Vg/U ~ —1,0, but a suppression of 
charge fluctuations in between. While the spin fluctua- 
tions at the particle-hole symmetric point increase with 
increasing U /F^ the charge fluctuations decrease at the 
same time. Still, the conductance is always perfect, since 
the scattering of conduction electrons screens the mag- 
netic moment, giving rise to the large spin fluctuations. 
We characterize charge and spin fluctuations via: 

{5Sl,,f = {{Sl,,f)~{Sl,,f- (17) 
(5n,)2 - (n^) - (n,)^ (18) 
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Fig. 6. One dot, t' /t = 0.4, U/t = A, Vg ^ -U/2. (a) Spin- 
spin correlations vs. site, (b) integrated spin-spin correlations 
5z(r) [see Eq. (20)] vs. distance r. DMRG data are shown for 
N = 31 (circles, olo), N = 32 (squares, ole), N = 33 (stars, 
ele); M — 300. See the text in Sec. 3.2 for a discussion. 



Properties (ii) and (iii) of the single- impurity Anderson 
model can nicely be seen in Figs. 5(a) and (b): olo and ole 
clusters exhibit a broad maximum in the spin fluctuations 
around the particle-hole symmetric point Vg = —U/2, 
while charge fluctuations are suppressed [Figs. 5(a) and 
(b)]. The maxima in (fcdot)^ vs. Vg/U at Vg/U = -1,0 
are due to the charge degeneracy between states with 
ndot = and 1 and rtdot = 1 and 2, respectively. For 
comparison. DMRG results for N = 200 sites (triangles) 
are included in Figs. 5(a) and (b). 

As for the ele clusters, we emphasize that the main 
consequence of the finite spin projection (S^^^) present for 
~1 ^ ^0 I^ig- 5(c)] on the dot is a reduction 

of spin fluctuations [see Fig. 5(b)]. A Kondo-resonance, 
however, which is at the heart of perfect transmission 
through a quantum dot at the particle-hole symmetric 
point, emerges as a consequence of virtual charge fluc- 
tuations due to scattering of conduction electrons off the 
impurity that cause substantial spin fluctuations on the 
dot, thus screening its moment [61]. 

Apart from spin fluctuations, and equally importantly, 
spin-spin correlations between the dot and the leads are 
affected as well. Figure 5(d) shows the spin-spin correla- 
tions between the dot and the first site in the leads: 



Sf. 



(19) 



which are substantially smaller around the particle-hole 
symmetric point Vg = —U/2 on ele clusters than for other 
clusters. These two features illustrate the failure of ele 
clusters to harbor precursors of Kondo physics as long 
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Fig. 7. One dot, t' /t = 0.4, U/t = 4, Vg = -U/2. Finite-size 
scaling analysis of static properties: (a) charge fluctuations on 
the dot; (b) spin fluctuations on the dot; (c) local spin (^^ot) on 
the dot (compare Ref. [62]); (d) spin-spin correlations between 
the dot and its neighboring sites. In (d), the upper ole curve 
is for i in the right lead (i — dot + 1), while the lower one is for 
i in the left lead {i = dot - 1). DMRG results {M < 500) are 
shown for olo (circles), ole (squares), and ele clusters (stars). 



as system sizes are small. We shall discuss the finite-size 
scaling below. 

Let us next discuss how the emergence of screening is 
reflected in the spatial extension of spin-spin correlations. 
DMRG results for spin-spin correlations [see Eq. (19)], 
measured away from the dot, and for their integral Sz{r), 



E 



5, 



r'{i) 



Mot 



(20) 



are displayed in Fig. 6(a) and (b), respectively, for = 31 
(olo), N = 32 (ole), and A^ = 33 (ele). We first observe 
from Fig. 6(a) that, while the leads are symmetrically po- 
larized in the case ofA^ = 31, ole clusters mainly polarize 
one lead, namely the one that has an odd number of tight- 
binding sites. In this case, an approximate singlet is formed 
with this lead. When integrated over distance, the spin 
correlations behave quite similarly for these two cluster 
types, i.e., ole and olo. In contrast to that and consistent 
with the notion of the absence of Kondo physics on small 
ele clusters, only small spin-spin correlations spread out 
into the leads for these clusters. 



3.3 Static properties: Finite-size scaling 

We proceed by discussing the finite-size scaling of static 
properties. Our DMRG results for system sizes 15 < A^ < 
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400 are displayed in Fig. 7. The main observations are: (i) 
spin and charge fluctuations converge the fastest on ole 
clusters; (ii) in the case of ele, the dot carries a large, finite 
(S^Qj) that decays slowly with TV, suppressing spin fluc- 
tuations and spin-spin correlations with the dot's neigh- 
boring sites; (iii) antiferromagnetic spin-spin correlations 
•^fdot ^'"^ enhanced in the case of olo clusters on small 
systems, while on ole clusters, mainly the lead with an 
odd number of sites develops substantial AFM correla- 
tions with the dot. 



3.4 ECA results for the conductance 

Early studies [27] using this technique have phenomeno- 
logically proposed that for one dot, the best results are 
obtained when it is neighbored by an odd number of sites 
(olo), while ele clusters exhibit a Coulomb-blockade like 
behavior. Using insights of Sees. 3.2 and 3.3, we can now 
explain this observation on a more rigorous basis. Small 
ele clusters typically locate the S^^^^^ = 1/2 to a large 
extent onto the quantum dot, thus giving rise to a split 
resonance. Consistently, we find a charge gap in the local 
density of states (LDOS). While the embedding procedure 
manages to redistribute spectral weight in the LDOS quite 
well [84] , it does not close the charge gap given by the ex- 
actly solved interacting region. 

For illustration, we show ECA results for one dot at 
t'/t = 0.4, U/t = 4 and aU three cluster types in Fig. 8(a), 
(c), and (d). The conductance is computed from Eq. (14). 
The figure also contains the conductance as obtained from 
the Friedel sum rule (FSR) using the charge computed 
with static DMRG [61]: 

G = Gosin2(7r(ndot)/2), (21) 

where Go ~ 2 due to spin degeneracy. Here we claim that 
the charge (ridot) can be obtained to high precision from 
ground-state DMRG. To support this we show a compari- 
son for the conductance derived from Eq. (21) with exact 
Bethe-ansatz resuhs [60] in Fig. 8(b) for U/T = 25. Note 
that in the wide-band limit U < At, the hybridization pa- 
rameter is = 27rt'^ Plead, and since we work with semi- 
infinite-leads. Plead = l/(7rt) [79]. Our DMRG results are 
for U = 0.5t, t'/t = 0.1. We find excellent agreement with 
the results from Ref. [60] in the limit of a large U/F, which 
justifies the use of static DMRG to benchmark the ECA 
results at smaller values of J7/I^ ^ 12. 

Clearly, the ole clusters are closest to the FSR as 
shown in the main panel of Fig. 8(a). In Fig. 8(d), we 
further observe a Coulomb blockade type-of behavior for 
ele clusters. Regarding the olo cluster, although it quali- 
tatively provides the correct physical description, the con- 
ductance plateau comes out too narrow as the charge 
around the particle-hole symmetric point Vg = —U/2 varies 
too fast as compared to the DMRG result [see Fig. 8(c)]. 
Note that the ECA results for the charge (ridot) as dis- 
played in Fig. 8 are obtained after embedding by integrat- 
ing over the imaginary part of the dot's onsite Green's 
function. 




1 -1 
V /U V /U 



g g 

Fig. 8. One dot, t'/t = 0.4, U/t = 4. (a) Conductance G and 
charge (jidot) vs. gate potential Vg/U for ole (ECA, Neu ~ 12, 
solid lines), and DMRG results for the charge and G obtained 
through the FSR (circles and squares, respectively). ECA re- 
sults for an _Zo cluster with A^ed = 10 sites are indistinguishable 
from the ole curve, (b): Comparison of conductance obtained 
from the Friedel sum rule using DMRG data for the charge 
(^^dot) (A = 300, M = 500, U/t = 0.5, t'/t = 0.1) with Bethe- 
ansatz results from Ref. [60]. (c) Conductance G and charge 
('^dot) vs. gate potential Vg/U for olo (ECA, Aed = H, solid 
lines), and DMRG results for the charge (circles), (d) Con- 
ductance G and charge (ridot) vs. gate potential Vg/U for ele 
(ECA, Aed = 9, solid lines), and DMRG results for the charge 
(circles). 



A further improvement of results can be achieved by 
using a transformation onto bonding and anti-bonding or- 
bitals in the Hamiltonian for one dot, Eqs. (l)-(3), and 
(15) (see, e.g., Rcfs. [10,11] and references therein): 

- ^[cL<. +4,.]; A,a = ^[4,. - • (22) 

Here, the index i measures the distance away from the 
dot and iL,iR are the corresponding indices in the two 
semi-infinite leads. The transformation decouples the dot 
from the anti-bonding states, leaving it coupled only to the 
bonding states. This is technically a great simplification as 
leads almost twice as large can now be treated. While such 
transformation cannot be used for the spatially asymmet- 
ric ole clusters, it works for ele and olo. We observe that 
exactly the same result is found for ele, comparing, say, 
the full A^ = 9 =4-dot-4 with the symmetrized dot-4 clus- 
ters. Here, the spin still cannot be screened in a singlet 
state. For olo, the results substantially improve, as now, 
effectively a dot-odd (lo) cluster is solved that has an over- 
all -S'^qj^ji = ground-state, allowing for singlet formation. 
This gives results similar to those obtained from an ole 
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configuration witli less computational effort, as a smaller 
cluster is diagonalized. 

Our observations lend strong support to the idea that 
when the global quantum numbers of the exactly solved 
cluster allow for a singlet state to be the ground state - 
as is the case for ole and lo clusters - then a precursor of 
the Kondo effect can be seen even on small ED clusters, 
with a favorably fast convergence with system size. 

Finally, Figs. 8(a), (c), and (d) suggest a way to gauge 
the quality of EGA data even in a situation where no ex- 
act results for the conductance are available. We sec that 
the cluster type that produces the best charge vs gate 
potential curve yields the best G vs. Vg curve as well. 
Thus, performing a finite-size scaling analysis of proper- 
ties such as the charge or spin fluctuations using either a 
Lanczos routine (which is a part of standard EGA codes, 
see Sec. 2.2) or DMRG is well suited to identify the opti- 
mum cluster type to be used in conductance calculations 
in EGA. 



3.5 One quantum dot: Summary 

Let us summarize our observations from this Sec. 3. We 
can understand the behavior of reflection-symmetric clus- 
ters with an odd number of sites and an S'^ta,! = 1/2 
ground state in terms of rigid spins (5'f ) pinned to certain 
sites. In the case of eie, the dot is mostly polarized and 
carries a finite magnetization. This is qualitatively equiva- 
lent to the application of a magnetic field that causes Zee- 
man splitting of the Kondo resonance and, thus, a suppres- 
sion of the conductance [61]. For the case of olo clusters, 
the leads are polarized, with only every second site car- 
rying a finite magnetization [62]. For comparable system 
sizes, a substantially larger current is observed on these 
clusters, which, in tDMRG simulations, however, does not 
exhibit a current that is constant in time for t' ^ t. 

Our analysis further suggests that either ole or lo 
clusters are best suited to study the emergence of Kondo 
physics with techniques such as tDMRG or EGA, while 
on system sizes accessible to these techniques, ele clus- 
ters yield qualitatively different results with essentially a 
Goulomb-blockade like behavior. The overall conclusion is 
that in order to provide a reliable description of the many- 
body effects at the Fermi level, clusters with a S^^^^i = 
ground state are necessary in calculations with these two 
methods. Among these, the advantage of io clusters is 
that larger system sizes can be accessed. Finally, we stress 
that no matter what the cluster type is, if one was able 
to go to extremely long chains, the same physics would be 
recovered [62]. 



4 Three Dots 

We now turn to a more complicated case, three serially 
coupled dots, as shown in Fig. 1(b). Apart from the in- 
terest in the emergence of Kondo physics, many quantum 
dots arranged in an array can be considered an interpola- 
tion between a single localized impurity and a bulk Mott 
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Fig. 9. Three dots, i! jt = 0.3, U jt = 1 (a) and U jt = 3 (b). 
(a) Spin projection (Sf ) vs. t" jt for the outer dots (squares) 
and tlie central dot (stars) for = 13 (o5o, dashed lines) and 
= 15 (e5e, solid lines) clusters, (b) Same as in (a), but for 
U jt = 3. Note that the magnitude of the fixed spins for both 
clusters (o5o and e5e) increase with the value of U . On o3e 
clusters, (Sf ) = (DMRG; M = 300). 



insulator. In the first case, the transmission at half filling 
is perfect, but zero in the later due to the presence of the 
Mott gap [85]. 

The interacting region is described by 



iJint = ~t" + ^■'^■) 

i=l 

3 3 

+uY^n,^^ni^^ + Y^Vgn,. (23) 



To simplify the notation, we will refer to the dots as Dl, 
D2, and D3, where D2 is the central dot. We focus on 
the behavior at the particle-hole symmetric point unless 
stated otherwise. In the case of three dots, we distinguish 
between odd-3-odd (o^o), odd-3-evcn (o^e), and even-3 
even (e5e) clusters. 

Following Refs. [40,67], for large t" > U/^/2, the three 
dots can be viewed as a molecule with, at half filling, 
two electrons occupying its lowest state, and a result- 
ing S't^otai = 1/2 from the third electron. For intermediate 
t" ^ 0.3J7, it has been suggested that the three dots be- 
have as a linear antiferromagnet due to strong AFM spin- 
spin correlations between them. Finally, at small t" ^ [/, 
the system enters a very subtle two-stage Kondo regime. 
First, and below a Kondo scale the outer dots form 
a Kondo singlet with their adjacent leads. Then, at tem- 
peratures below a second Kondo scale T < T^^ < 
the central dot is screened by the quasi-particles of the two 
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Fermi liquids formed by the outer dots and their leads. Al- 



though 



(2) 



is orders of magnitude smaller than T, 



(1) 

K 



"40, 



86] , which renders an experimental confirmation of the ef- 
fect a challenging task, nevertheless the two-stage process 
conceptually plays an important role for extremely low 
temperature physics. In all three cases and at tempera- 
ture T = 0, the system possesses a spin 5 = 1/2 Kondo 
ground state. 

Concerning the issue of the conductance, several stud- 
ies using a variety of approaches such as perturbation 
theory [65], NRG [67,68], and functional renormalization 
group [36] report perfect conductance at the particle-hole 
symmetric point due to the Kondo effect. An earlier EGA 
work on this model [28] is at odds with this conclusion. 
In particular, an exact zero of the conductance at the 
particle-hole symmetric point has been found for all values 
oit" studied, i.e., t" < 0.5t. 

We show that the vanishing of G at V^, = —U/2 at 
t" /t < 0.3 is due to the choice of o3o clusters for the EGA 
calculations of Ref. [28]. At large t" > t', this dip is an 
artefact of certain approximations taken in the embedding 
procedure: if we perform the EGA according to the pre- 
scription of Sec. 2.2 without mixing Green functions for 
spin up and and down before the embedding is carried 
out, this artifact will disappear on o3o clusters. 

We shall outline below that the optimum cluster type 
- in the sense of fastest convergence of static properties 
with system size - is o3e. On this cluster type, no spuri- 
ous dip at Vg = —U/2 exists, independent of the spin- 
mixing strategy adopted in EGA calculations, and the 
EGA results then are in agreement with the conclusions 
of Refs. [36,65,67,68]. 

The section is organized as follows: We first analyze 
static properties and their finite-size scaling in Sec. 4.1. 
Section 4.2 contains our EGA results, obtained following 
the procedure discussed in Sec. 2.2, for the conductance. 
In Sec. 4.3, we briefly discuss the results one obtains in 
the regime of t" < t' using o3e clusters. 



4.1 Static properties: Finite-size scaling 

Let us first consider the three dots isolated from the leads 
(t' = 0) at [/ = i as a function of t" /t, which will help 
to understand e3e and o3o clusters. In the S^^^^i = 1/2 
subspace, the two lowest-lying states are both doublets 
{S = 1/2), however, with different spin projections (Sf) 
on the dots. The ground state has - schematically - an 
up-down-up pattern, while the first excited state realizes 
an up- up- up structure, with most of the ^totai = 1/2 on 
the central dot. The actual distribution of (5,?) depends 
on t" /U for a given t', but when the isolated cluster is em- 
bedded in leads with an overall odd number of sites (i.e., 
all sites in the leads plus the three dots), the two low- lying 
states are mixed with different weights: on o3o clusters, 
an down-up-down configuration is preferred, while on e3e 
ones, the original ground-state pattern (up-down-up) pre- 
vails. We illustrate this in Fig. 9 where we plot (Sf) vs. 
t"/t' for = 13 {o3o) and iV = 15 (ele). We further see 
that, as a function of increasing t" /t, the spin projection 
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Fig. 10. Three dots, t' /t = 0.3, U/t = 1. Finite-size scaling 
of spin-fluctuations on the dots for t" /t — 0.05 [(a),(b)] and 
t"/t = 1 [(c), (d)]. DMRG resuhs for aU three cluster types: 
o3e (squares), eSe (circles), and o3o (triangles) [M < 500]. 



on the central dot (displayed with stars in the figure) is 
substantially reduced and moved to the leads in the case 
of o3o. The large (Sf ) located on the outer dots (displayed 
with squares in the figure) seen in the case of eSe decays 
much slower as a function of t" /t. 

The rigid spins on either of these cluster types imply 
suppressed spin fluctuations for the individual dot (D2 in 
the case of o3o and Dl, D3 in the case of e3e) as well 
as modified spin-spin correlations. As two case examples, 
we present the finite-size scaling of such quantities with 
system size in Figs. 10 and 11 for U ~ t, t' /t ~ 0.3 with 
t"/t = 0.05 and t"/t = 1. 

Let us first discuss e3e clusters (circles in Figs. 10 and 
11). Due to the large spin projection on the outer two dots, 
spin fluctuations as well as the spin-spin correlations with 
the first site in the leads are suppressed. The latter is an 
indicator of how strongly the dots polarize the leads. Since 
the spin-projections Sf^^ and Sf-,^ are - both at small and 
large t" - only slowly redistributed to other sites as the 
system size grows, as is evident from Figs. 10(a) and (c), 
we conclude that these clusters will produce low conduc- 
tance values or currents in EGA and tDMRG calculations, 
respectively. 

The o3o clusters (triangles in Figs. 10 and 11) are the 
ones originally used in Ref. [28]. Figure 10(b) reveals that 
spin fiuctuations on the central dot are suppressed in the 
limit of small t" < t, which is due to a large fraction 
of S^^^^^ located on the central dot. This cluster type is 
therefore geometrically similar to the ele clusters of one 
quantum dot. As it is illustrated for the case of t"/t =- 
0.05 in Fig. 10(b), the spin fiuctuations on the central dot 
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of spin-spin correlations for t" /t = 0.05 [(a),(b)] and t" /t = 1 
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tions between the central dot and the outer ones. In the o3e 
case (squares), there are two branches as in reality, we alter- 
nate between o3e and eSo (even-dots-odd) configurations as A'' 
increases, (b) and (d): spin-spin correlations between the outer 
dots and the first site in the adjacent lead. In these two panels 
and for o3e clusters (squares), the upper(lower) curve shows 
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increase very slowly with system size. We may therefore 
argue that here the central dot fails to participate in a 
Kondo effect, especially in the two-stage Kondo regime. 
This pathology, i.e., the finite S^2^ the origin of the 
conductance dip at Vg = —U/2 observed with EGA in 
Ref. [28], similar to the case of one quantum dot. 

Clusters with an overall S^^^^i = (i.e., o3e - squares 
in Figs. 10 and 11) exhibit the least significant finite-size 
effects and local fluctuations converge the fastest with sys- 
tem size N, while spin-spin correlations exhibit an asym- 
metry between the left and right lead. As a conclusion, 
we may expect reliable results for conductances from this 
cluster type, and this notion will further be corroborated 
by EGA results in Sec. 4.2. 



(a) U=t, t7t=0.3, t"=t; ECA: o3e 




1.2 



0.8 



1.2 

a°o.8 
a 

0.4 


1.2 

a° 0.8 
o 

0.4 


1.2 

0° 0.8 

0.4 





Fig. 12. Three dots, t' /t = 0.3, U/t = 1, t" /t = 1. Con- 
ductance and charge vs. gate potential Vg/U. (a) ECA results 
for the conductance and charge for an o3e cluster (solid lines, 
A^ED = 12) vs. DMRG results for the charge ntotai (dotted line, 
N = 300, M = 400) and NRG results ([87], circles), (b) ECA 
results for an o3o cluster (solid lines, Ned = 9). (c) ECA re- 
sults for an e3e cluster (solid lines, Aed = 11)- DMRG results 
for ntotai (dotted line) are included in (b) and (c) for compar- 



and (c) (eSe). As expected from the discussion of static 
properties in Sec. 4.1, the results for ritotai for o3e clusters 
obtained after embedding agree very well with DMRG re- 
sults [squares, see Fig. 12(a)], while the results for the con- 
ductance from other clusters deviate strongly from DMRG 
data [squares, see Figs. 12(b) and (c)]. From Fig. 12(c) 
[solid lines], we realize that eSe clusters produce the ex- 
pected Goulomb-blockade like behavior due to the rigid 
spins on the outer dots (see Fig. 9). 

Figure 12(a) also contains NRG data from Ref. [87]. 
The agreement between EGA and NRG [87] is quite good, 
especially for the central Kondo peak and the Goulomb 
blockade valleys at ntotai = 2,4. This clearly estabhshes 
that reliable results can be obtained from EGA calcula- 
tions using o3e clusters. 



4.2 ECA results for the conductance of three dots: 
Comparison of cluster types 

4.2.1 Gate-potential dependence 

In this section, we present EGA results for the conduc- 
tance G and ritotai ("total = (n-Di) + ("02) + ("Ds)) for the 
three dots configuration as a function of the gate poten- 
tial, comparing all three cluster types. The EGA results 
(solid fines) are depicted in Fig. 12(a) {o3e), (b) (o5o). 



4.2.2 Behavior at the particle-hole symmetric point 

Thus far, we have seen that with o3e clusters, we obtain 
the best agreement with NRG regarding the gate-potential 
dependence. Let us next address the value of the conduc- 
tance at the particle-hole symmetric point as a function 
of t" . This will indicate what parameter range can be ac- 
cessed when decreasing t" /t using o3e clusters. 

Figure 13 compares o3e and o3o EGA conductance re- 
sults at Vg = —U/2 for U = t and U = 3t. In the former 
case, where o3o clusters with A'ed = 5,9 arc compared 
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Fig. 13. Three dots with t' /t = 0.3, U/t = 1, and U/t = 3; 
Variation of the conductance computed with ECA with t" at 
the particle-hole symmetric point Vg = — [7/2. The results for 
U/t = 1.0 are calculated using o3o with N — 5 sites (circles) 
and A'^ = 9 sites (triangles), as well as o3e clusters with — 
12 sites (squares). At this value of U/t, there is barely any 
difference between the two types of clusters, i.e., as the value 
of t" increases, the conductance at the particle-hole symmetric 
point changes from to Go at basically the same value of t" , 
independent of cluster type. However, the same is not true for 
the U/t = 3 results (iV = 10 for o3e cluster (stars) and N = 9 
for o3o cluster (diamonds)). It is clear that the o3o clusters 
yield a vanishing G at substantially higher t" than the o3e 
ones. As explained in the text, this is caused by a larger fixed 
spin in o3o clusters at higher U values. 



with an oSe cluster with A^ed = 12, it is shown that there 
is barely any difference between the results for the two 
different types of clusters: for t" /t « 0.5, the correct re- 
sult G/Go = 1 is recovered. However, in the later case 
{U = 3t), there is a pronounced difference between o3o 
{Ned = 9) and o3e {Ned = 10) results: the vanishing 
of the conductance at the particle-hole symmetric point 
for the o3o cluster occurs at a much higher value of t" /t 
than for the o3e cluster. The reason for that can be under- 
stood by analyzing the results in Fig. 9, where the fixed 
net-spin present in o3o clusters is compared for U = t 
[panel (a)] and U — 3t [panel (b)]. A larger value of U, 
for a fixed t" , increases the magnitude of the fixed spins 
in the three dots (see Fig. 9). Therefore, as pointed out 
above, o3e clusters, with S^^^^^ ~ 0, provide better results 
at the particle-hole symmetric point when compared to 
o3o clusters. In the next subsection, we will discuss why 
the conductance eventually also vanishes for low t" values 
in the o3e cluster and numerically accessible values of L, 
which do not exhibit finite values of (S^). 

We now turn to a finite-size scaling analysis of G{N) 
at Vg = —U/2. Using cluster sizes within ECA that are 
numerically tractable (i.e., A^ed ^ 13), we find that, at the 
particle-hole symmetric point, for U/t = 1.0, G(A^ed)/Go 
scales to one with I/A^ed for t"/t > 0.3 (upper panel in 
Fig. 14). In addition, it is clear from the results shown in 
Fig. 14 (b) that, the larger the value of U, the higher is the 
value of t" below which one cannot clearly ascertain the 
extrapolation of G{Ned) to Go at small t" values. Obvi- 
ously, this happens because the Kondo temperature of the 
second stage decreases with U, causing the ECA conduc- 
tance to vanish at a larger t" . Summarizing, both o3o and 



o3e clusters produce results in qualitative agreement with 
other techniques at Vg ~ —U/2 above certain values of t". 
For U/t = 1, G = Go can be obtained by extrapolation 
with a G{Ned) = Go - const /A^ed function for t"/t > 0.2. 
For smaller t", system sizes are too small to establish a 
clear trend of G{Ned)/Go towards unity due to a finite 
level spacing. 

We are led to conclude that ECA captures quite well 
the correct physics depending upon the parameters and 
the cluster taken. The optimum cluster can be identified 
by analyzing the finite-size scaling of, e.g., the charge, with 
either the Lanczos solver that is part of ECA or indepen- 
dent techniques such as ground state DMRG, the latter 
being available in open access releases, of, e.g., the ALPS 
project [88]. 



4.3 o3e clusters in the small t" regime 

We next wish to illustrate the finite-size effects encoun- 
tered in the small t" regime, for o3e clusters. Unfortu- 
nately, when t" is decreased, the asymmetry between the 
two outer dots on o3e becomes quite pronounced. Note 
that the total charge from ECA still agrees well with 
DMRG down to very small t" (not shown). However, spin- 
spin correlations on small o3e clusters suffer from (i) the 
aforementioned asymmetry between the two outer dots 
(see Fig. 11), and (ii) the fact that the central dot is ba- 
sically decoupled from the leads (at the Fermi energy). 
This is indicated by small spin-spin correlations between 
the central dot and the first site in the leads, which indi- 
cates that the effective exchange with the leads is smaller 
than the energy level spacing. Consistent with the first 
point discussed here, in accordance to the LDOS calcu- 
lated with ECA for an o3e configuration (not shown), 
at small t" and for Vg = ~U/2 (particle-hole symmet- 
ric point), only one of the outer dots develops a sharp 
Kondo resonance, namely the one connected to the odd 
lead. The dot connected to the even lead can only be 
Kondo screened through the other two dots. However, at 
low values of i", this mechanism becomes ineffective and 
thus the dot's LDOS at the Fermi energy is low, leading 
to the suppression of the conductance at the particle-hole 
symmetric point, as seen in Fig. 13. Therefore, even for 
the optimum cluster type o5e, the two-stage Kondo regime 
seems to be out of reach within the real-space variant of 
ECA and current computational resources. 



5 Summary 

In this work, we have performed an extensive finite-size 
scaling analysis of fluctuations and correlations in nanos- 
tructures, such as one and three quantum dots coupled 
to noninteracting, tight-binding leads. Strong differences 
in the finite-size scaling behavior emerge depending on 
even-odd effects [62]. These findings are relevant for the 
interpretation of results from numerical approaches for 
the calculation of the conductance of nanostructures that 
are based on analyzing clusters of finite length with open 
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Fig. 14. Three dots, t' /t = 0.3, U/t = 1 (a) and U/t = 3 (b). 
(a) Finite-size scaling analysis of ECA results on o3e clusters 
for G{Neu) at the particle-hole symmetric point {Vg — —17/2) 
for t" /t = 0.45,0.4,0.3,0.2 (squares, triangles, circles, and 
stars). Dashed lines are the linear fits to G{Neu)- (b) Same 
as in (a), but now for U/t = 3, and t" /t = 0.8, 0.7, 0.6, 0.5, and 
0.4 (circles, triangles, diamonds, stars, and squares). Note that 
in this case, the threshold, below which the convergence to Go 
cannot be ascertained, is higher than in the U/t = 1 case. 



boundary conditions, such as time-dependent DMRG or 
the embedded-cluster approximation (ECA). To qualita- 
tively capture precursors of the Kondo effect within these 



two methods, clusters with a S^^^^l = are best suited, 
as electrons on the dot at half filling naturally partici- 
pate in a singlet state. Other cluster types with an over- 
all S, 



total 



1/2 may be far away from capturing Kondo 
physics in the sense of slow convergence with system size. 
For certain configurations, most of the S^^^^^ = 1/2 is 
found on the central site, which, if this site is a quan- 
tum dot, will cause a gap in the local density of states at 
the Fermi level, reminiscent of the Zeeman-splitting of the 
conductance due to the application of a magnetic field. 

These findings help to understand the discrepancy be- 
tween ECA results for three dots and those of other tech- 
niques: the vanishing of the conductance at the particle- 
hole symmetric point obtained with ECA is only seen on 
certain cluster types, but disappears when the three dots 
are embedded in a cluster with an overall S^^^^^^ = 0. In 
the latter case, the ECA results are then in agreement 
with the picture proposed in Refs. [36,40,65-69]. While 
parameter regimes with exponentially small energy scales 
such as two-stage Kondo regimes are difficult to be ac- 
cessed with this technique in its real-space variant, very 
good agreement with exact results such as the Friedel sum 
rule in the case of one dot - as shown in Fig. 8(a) - or 
NRG in the case of three dots is otherwise found in other 



regions of the parameter space. The latter is established 
from Fig. 12(a) and Fig. 13. 

Another model that exhibits a two-stage Kondo ef- 
fect is two dots coupled in a T-shape geometry [86,89] At 
half filling, the conductance is found to vanish [86,89], a 
picture that disagrees with ECA results in the two-stage 
Kondo regime [79]. We have performed a similar analy- 
sis of static properties and ECA results on different clus- 
ter types for this geometry, with results not shown here. 
While in the regime of large t" /U, where t" is the hop- 
ping between the two dots, good agreement between the 
Friedel sum rule [86] and ECA is found, the two-stage 
Kondo regime, which emerges at small t" [86,89] suffers 
from strong finite-size effects. Unfortunately, no cluster 
type shows a particularly fast scaling behavior here. 

We conclude that the determination of the optimum 
cluster type is crucial. This can be done by studying dif- 
ferent cluster types with the Lanczos solver that is part 
of standard ECA codes. Further, a comparison of ECA 
results for static properties, such as the total charge of a 
nanostructure as a function of gate potential and very im- 
portantly, spin and charge fluctuations, against indepen- 
dent techniques, such as ground-state DMRG, is helpful 
in this respect as well, as we showed in this work. Promis- 
ing results from a new version of ECA that incorporates a 
logarithmic discretization of the DOS of the leads clearly 
indicate that the range of validity of the embedded-cluster 
approximation can further be improved [84] . 
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